LA-UR-12-24539 



Precise determination of the structure factor and contact in a unitary Fermi gas 

Sascha Hoinka, 1 Marcus Lingham, 1 Kristian Fenech, 1 Hui Hu, 1 

Chris J. Vale, 1 Joaquin E. Drut, 2,3 and Stefano Gandolfi 3 

Centre for Atom Optics and Ultrafast Spectroscopy, 
Swinburne University of Technology, Melbourne 3122, Australia 
2 Department of Physics and Astronomy, University of North Carolina, Chapel Hill, NC 27599-3255, USA 
3 Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

(Dated: September 19, 2012) 

We present a high-precision determination of the universal contact parameter in strongly inter- 
acting Fermi gases. The contact governs the high-momentum properties (i.e. the short-range corre- 
lations) of these systems. The experimental measurement utilises Bragg spectroscopy to obtain the 
dynamic and static structure factors of ultracold Fermi gases at high momentum in the unitarity and 
molecular Bose-Einstein condensate (BEC) regimes. We have also performed quantum Monte Carlo 
calculations of the static properties, extending from the weakly coupled Bardeen-Cooper-Schrieffer 
(BCS) regime to the strongly coupled BEC case, which show excellent agreement with experiment at 
the level of a few percent. At unitarity our low temperature measurement provides a new benchmark 
for the zero temperature homogeneous contact which we find to be 3.33 ± 0.07. 
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Ultracold atomic gases are unparalleled as a means to 
quantitatively probe strongly coupled systems which lie 
at the intersection of condensed-matter pQ , nuclear [51 [3] 
and high-energy physics [HE]. In this context, two- 
component mixtures of fermionic atoms have particular 
significance as they are generally stable against inelas- 
tic decay. These universal quantum systems [3HS] are 
characterized by strong coupling in the form of s-wave 
interactions of short range r and large scattering length 
a, such that the only scales left are thermodynamical: 
the density n or chemical potential [i and the tempera- 
ture T, just as in a non-interacting gas. This situation, in 
particular the unitary limit <— k F r <1< k F a — > oo, 
where k F is the Fermi wavevector [3], represents a ma- 
jor theoretical challenge as there arc no small parameters. 
While there is no exact description, a variety of approx- 
imate techniques have been developed; however, these 
often predict quite different outcomes. 

One of the key quantities characterizing these sys- 
tems is the universal contact parameter C, introduced by 
Tan [TUl [TT] . The contact derives from the short-range 
properties of strongly interacting quantum gases and is 
the cornerstone of a number of exact relations describing 
the static and dynamic properties such as the equation of 
state and dynamic response functions [12-15 . Evaluat- 
ing these exact relations requires precise knowledge of C 
itself, which is very challenging to compute with different 
calculations varying by as much as 10% (TH [T7] . 

In this letter we provide a new benchmark measure- 
ment, accurate to the level of 2%, for the contact at 
unitarity. This measurement is furnished by a precise 
determination of the dynamic and static structure fac- 
tors using Bragg spectroscopy. In addition, we present 
new Quantum Monte Carlo (QMC) calculations accurate 
to the level of a few percent. Our results indicate that 



theory and experiment are reaching a new level of con- 
vergence, showing that this difficult problem is becoming 
tractable. 

Experimental setup. The experiments presented here 
use a gas of 6 Li prepared in an equal mixture of the \F = 
1/2, rriF = ±1/2) spin states, evaporatively cooled in a 
single-beam optical dipole trap. Interactions are tuned 
to the unitary limit by setting the magnetic field to 833.0 
G, close to the pole of the Feshbach resonance. Following 
evaporation, the cloud of approximately N/2 = 2.2 x 10 5 
atoms per spin state is at a temperature of ~ 0.07Tp 
where Tp is the local Fermi temperature in the centre of 
the trap. The atoms are then loaded into a second optical 
dipole trap produced by a 10 W single frequency 1064nm 
fibre laser that is spatially filtered before entering the cell 
to produce a deep trap with large harmonic region. The 
trapping frequencies in the radial and axial directions are 
w x = uj y = 2nx97 Hz and lu z = 27rx24.5 Hz, respectively. 

Bragg scattering is performed as in previous work |18l 
I19j . Two laser beams, detuned by approximately 600 
MHz from the nearest atomic transition, illuminate the 
atom cloud, intersecting at an angle of 26 — 84°. This 
sets the probe wavevector k = (Aw /X) sin 8 where A = 
671nm. The beams are obtained from the same laser 
with two different acousto-optic modulators, driven by an 
amplified signal from a multichannel direct digital syn- 
thesizer. This allows us to precisely tune the frequency 
difference between the two beams u) with a stability bet- 
ter than 1 Hz. 

To obtain the dynamic response function we measure 
the momentum imparted to the cloud for a range of Bragg 
frequencies. By imaging atoms in both spin states with a 
short (850 /is) time delay between images we can deter- 
mine the centre of mass cloud displacement AX, which 
is proportional to the imparted momentum |18j , without 
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FIG. 1. (color online) Density response functions of a strongly 
interacting Fermi gas (a) in the unitary limit (833 G) and 
(b) for fcp°a = 1 (783 G), as a function of lu/lj,,. These are 
equivalent to the dynamic structure factor S(k, ui) at the high 
momentum used for the measurement (k — 4.46fcp°). 



being sensitive to fluctuations in the initial cloud posi- 
tion. We also perform a careful check of the response of 
the cloud at various Bragg frequencies in order to deter- 
mine the maximum Bragg intensity we can use and re- 
main in the linear response regime |19j . We can therefore 
maximize the signal-to-noise ratio at any given Bragg fre- 
quency by using the highest beam intensity that is within 
the linear region. Scaling the data by the product of the 
Bragg beam intensities thereby allows us to combine data 
measured at different intensities into a single spectrum. 
With these techniques, and by averaging approximately 
ten measurements at each w, we dramatically improve 
our measurements of the dynamic response. 

Figure [T] shows Bragg spectra obtained (a) at unitarity 
and (b) at fc^a = 1 where kf° = [2m/iw(3A0 1/3 ] 1/2 is 
the Fermi wavevector in a harmonic trap, m is the 6 Li 
atomic mass and Co = (w^cl^u^) 1 / 3 . These spectra con- 
tain a narrow peak at w r /2 arising from the scattering 
of pairs and a broader single atom peak centred around 
uj r [TH]. Error bars are the statistical standard devia- 
tion of the data obtained at a particular frequency. Each 
spectrum is scaled by its first energy-weighted moment 
providing a normalized measure of the density response. 
At the high momentum used here (k = 4.46fc^°) the nor- 
malized response is equivalent to the dynamic structure 
factor S(k,u>) [H1[5U]. The normalization exploits the /- 
sum rule for S(k,uj), which means the quantity plotted 
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this normalized spectrum over u provides the static struc- 
ture factor [21] which we find to be S(k = 4A6kf°) = 
1.182 ± 0.004 at unitarity, and 1.48 ± 0.01 at kf°a = 1. 
This represents a dramatic improvement over previous 
determinations which had an error bar on the order of 
10% [H]. 

With this precise determination of S(k) we can also ob- 
tain a new measure of Tan's universal contact parameter 
I for a trapped gas at unitarity. At this high momentum, 
S(k) = 1 + S n (k) PH [50], such that the dimensionless 
contact can be found directly as [2"2"] 
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is the dynamic structure factor in units of ui r /N where 
Lu r = lik 2 /(2m) is the recoil frequency. The integral of 



Using this expression we find I/(Nk^°) = 3.24 ± 0.06 
at unitarity and 12.5 ±0.2 at k^°a = 1. This is far more 
accurate than previously published data [121 H2] and 
sets a new benchmark for theoretical calculations of this 
quantity. This value is consistent with the one obtained 
from measurements of the equation of state |23] and from 
the frequency of collective oscillations [24]; however, it is 
significantly higher than what was found in both radio- 
frequency spectroscopy [55] [55] and photo-association 
data [571 EH]. Comparing with different theoretical re- 
sults our measurements are consistent with a Nozieres- 
Schmitt-Rink (NSR) calculation [16] and slightly below 
the zero temperature QMC results [33] but well above 
earlier many-body t-matrix calculations [30, 31j . On the 
BEC side our measured contact 12.5 ± 0.2 is slightly be- 
low the theoretical predictions ~ 13.2 possibly due to 
higher cloud temperature. 

Quantum Monte Carlo calculations. We compute the 
ground-state properties of strongly interacting fermions 
in the thermodynamic limit by means of QMC tech- 
niques. We calculate for a system of 66 fermions, which 
has been shown to be large enough to correctly describe 
the thermodynamic limit [32 , 33 . The QMC method has 
also proven to be sufficiently accurate in computing the 
energy as well other properties related to the contact in 
the unitary limit |29j . In this work, we use the same tech- 
nique to calculate the dimensionless energy per particle 
£ = E/Epc (also known as the Bertsch parameter) as 
a function of (k F a)~ l , where EpQ is the corresponding 
energy for the non-interacting gas, k F = (37r 2 n) 1//3 and 
n is the density. For each value of (fcpcs) -1 we first per- 
form a variational optimization of the many-body wave 
function, which we choose to be of the same form as in 
Ref. [29l [34] . The best variational ansatz is then used as 
trial wave function in the projection in imaginary time. 
The fixed-node approximation is used to control the sign 
problem, and the accuracy of the energy and other prop- 
erties depends on the quality of the variational wave func- 
tion. The fixed-node energy at unitarity is within a few 
percent with respect to the exact calculation of Ref. [33] . 

For each coupling strength we perform a different QMC 
calculation, varying the effective range r e of the two-body 




FIG. 2. (color online) The trap-averaged contact parameter 
as a function of fcp°a. Green solid line is obtained from QMC 
data combined with the LDA and the black dotted line is a 
Nozieres-Schmitt-Rink (NSR) calculation [TB|. Red circles are 
the new experimental values and the grey squares are the ex- 
perimental data from Ref. 22 . The inset shows the equation 
of state: points are the QMC simulations and the solid line 
is a functional fit to these. At unitarity, the experimental re- 
sult is 1/Nkp° = 3.24 ± 0.06, to be compared with the LDA 
value (QMC-based): 3.336 and the NSR result: 3.26 [16]. At 
fcp°a = 1 the experiment yields 12.5 ± 0.2, while the LDA 
result is 13.17, and the NSR value is 13.20 [16]. 



interaction to extrapolate to the r e — > limit. This 
is necessary because the energy per particle £ can be 
strongly dependent on r e , depending on the value of k F a. 
For (fcpCi) -1 > 0.2, we find that the slope of E(r e ) is 
negative and drops quickly as (fc F a) _1 is increased. Ex- 
trapolating to the r e — > limit is therefore crucial in the 
BEC region, while it is much less important on the BCS 
side of the Feshbach esonance. 

The extrapolated value of £ as a function of (k F a) -1 is 
shown in the inset of Fig. [2] One may then use this equa- 
tion of state to calculate the contact for the homogeneous 
system: 
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Using the above equation we find a value for the contact 
of 3.39 at unitarity and 12.78 for l/k F a = 1. 
The static structure factor is computed via 
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The QMC results for the homogeneous S(k) are 
shown in the main panel of Fig. [3] for [kpo)^ 1 = 
1, 0.5, 0.2, 0, -0.2, -0.5, -0.75 (curves from the top to 
the bottom in the figure), calculated at a finite effective 
range r e kp = 0.056. We have calculated the spin-parallel 



FIG. 3. (color online) Homogeneous static structure factor 
S(k) as a function of k/k F , for various coupling strengths 
computed using QMC. In the inset we compare the trap- 
averaged calculation of S(k) at unitarity with the present ex- 
perimental result (red diamond). The measurement is just 
over 1% higher than the QMC calculation. 



component of S(k), and found that Sff = iSj.j, = 0.5 
within the error bars for k > Akp in agreement with [20] 
and the experimental measurement of [T5] . 

Local density approximation. Knowing the homoge- 
neous contact and structure factor we can compute the 
trap-averaged values for these quantities using the local 
density approximation (LDA). The overall chemical po- 
tential /x is related to the local chemical potential fj,(r) by 
fi = fi{r) + V(r), where V(r) is the trapping potential. 
At unitarity, 
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Here s F (r) is the local Fermi energy, n(r) is the density 
profile, and we have assumed that £ does not depend on 
r (valid in the unitary and free-gas limits); more gener- 
ally we have £ = £(fcp(r)a), where k F (r) = (3-7r 2 n(r)) 1 / 3 
and a is the scattering length. For a harmonic trapping 
potential V(r) — ^mui 2 r 2 , the density profile is given by 



ri(r) = n(0) 



1 - 



R 2 



3/2 



R= [ie F (0)- 



1/2 



(6) 

and n(0) is determined by the total number of parti- 
cles in the system through the normalization condition 
J d 3 r n(r) = N. Given £, N and the frequency of the 
trap u, we completely determine the ground-state den- 
sity profile. Within the LDA, the total contact in a trap 
is given by I = J d 3 r X(r), where 
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and C is obtained from QMC data using Eq. ^. Using 
the density profile above, the contact in the trap is 
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In Fig. [2] we show the contact obtained from QMC simu- 
lation (solid green line), compared with the new experi- 
mental data presented in this letter (red circles) and the 
experimental data of Ref. [22] (grey squares). The LDA 
for the static structure factor in a trap is 
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where the dimensionless function S(k/k F ) was deter- 
mined via QMC. The inset of Fig. [3] shows the trap- 
averaged QMC result at unitarity. Also shown is the 
experimental measurement which lies just over 1% above 
the theory. Note that the LDA is less precise for 
S(k/k^°) than for the energy (and contact) as fits to 
the QMC data are less reliable for S(k), and because 
S(k) is computed at a small but finite effective range. 
Away from unitarity £ = £,{k F a), and 
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which follows from E = —PV + fiN combined with E — 
|Py— 8 ^ ra C. As in the unitary case, the central density 
determines /x, and the LDA equation /i = /i(r) + V (r) de- 
termines the density profile, which is solved numerically 
as k F d depends on r through k F (r). 

Homogeneous zero-temperature contact. While our 
measurements were performed on a trapped (inhomo- 
geneous) cloud our result at unitarity also provides a 
strong constraint on the zero temperature homogeneous 
contact. At finite temperatures the equation of state 
for the unitary Fermi gas is not known exactly, but for 
T <C T c , where T c is the superfluid transition temper- 
ature, trap averaged measurements are only weakly af- 
fected by the small population in the high temperature 
wings. Comparing calculations of the finite temperature 
trapped contact based on different models for the equa- 
tion of state QI1II3GIT] one finds that the ratio 1(T)/X(0) 
has a very weak temperature dependence in the range 
< T/T F < 0.08 with a relative discrepancy of less than 
1% for the different models. Thus, with knowledge of the 
trapped contact at a temperature T(<C T c ), it is straight- 
forward to obtain the zero temperature trapped contact 
with relatively low error. 

This procedure relies on a precise determination of 
the temperature which we achieve by fitting the finite 
temperature equation of state for the pressure [35] to 
the measured pressure obtained from one-dimensional 
(doubly integrated) density profiles [3H1 [37J . With this 
method we find a temperature of T/Tp = 0.067 ± 0.010. 
We then use the zero temperature relation Eq. (8) which 



connects the trapped contact to the homogeneous contact 
Co [H] and find cq = 3.33 ± 0.07 for the zero tempera- 
ture homogeneous contact density, in a very good agree- 
ment with the QMC result [25] . The increased error bar 
is primarily due to the uncertainty in the temperature 
measurement. We have used £ = 0.376 ± 0.004 [35J but 
note that the uncertainty in £ barely impacts the overall 
uncertainty as it appears in Eq. (8) to the 1/4-th power. 
An additional uncertainty arises due to the error bar on 
the precise magnetic field of the Feshbach resonance [33] • 
Based on a recent determination |40j . and our measure- 
ment performed at 833.0 G, our unitarity result for the 
contact would be shifted up by 0.07. 

In summary we have presented a high precision de- 
termination of the low temperature dynamic and static 
structure factors and contact of a strongly interacting 
Fermi gas. These systems are an ideal testbed for vali- 
dating different many-body calculations where exact pre- 
dictions are not available. Our measurements are now at 
a level that can discriminate between several of the es- 
tablished predictions and agreement with the latest QMC 
calculations is at the level of 1%. The measurement at 
unitarity also provides a new benchmark, accurate to the 
2% level, for the T — > limit of the homogeneous contact 
density which complements a recent measurement [2 6) at 
higher temperatures. 
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